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1  Introduction 


The  growth  of  films  on  single  crystal  surfaces  is  of  considerable  tech¬ 
nological  interest  which  has  been  discussed  extensively  in  the  literature  [1]. 
The  occurence  of  phase  transitions  in  the  adsorbed  layers  has  been  known 
for  some  time,  but  the  studies,  both  experimental  and  theoretical,  have  been 
done  in  either  dilute  gas-  solid  interfaces,  or  vacuum-  solid  interfaces  [2.  3], 
The  films  adsorbed  on  single  crystal  surfaces  do  exhibit  a  rich  variety  of  phase 
behavior[2].  There  is  good  reason  to  think  that  such  a  rich  phase  behavior 
also  exists  in  liquid  -  solid  interfaces.  However,  the  experiments  to  detect  it 
are  much  more  difficult  to  perform,  because  most  of  the  structural  surface 
analysis  is  carried  out  in  ultrahigh  vacuum,  and  the  observation  of  adsorp¬ 
tion  isotherms  is  possible  only  for  materials,  such  as  graphite,  that  have  a 
large  specific  surface,  and  therefore  very  irregular.  The  situation  in  electro¬ 
chemistry  is  different,  because  there  exists  a  technique,  polarography,  that 
is  surface  sensitive,  even  for  single  planar  surfaces,  and  which  is  capable  of 
detecting  the  sudden  formation  of  films  at  electrodes.  It  has  been  possible  to 
perform  structural  analysis  of  underpotential  deposits  of  metallic  monolayers 
[4]  on  metal  electrode  interfaces  using  X-ray  diffraction  techniques  [5,  6.  7], 
which  yield  structural  information  about  ordered  phases  of  adsorbed  mono- 
layers  on  single  crystal  surfaces.  The  underpotential  deposition  (UPD)  of 
monolayers  on  electrode  surfaces  occurs  suddenly  and  should  be  considered  as 
a  phase  transition  at  the  surface  [8,  9].  The  adsorption  isotherms  and  volta- 
mograms  of  systems  like  the  UPD  films  of  Cu  on  Au(lll),  Cu  on  Pt(lll). 
Pb  on  Ag(lll)  show  very  clearly  the  sudden  formation  of  ordered  films.  The 
crystal  nature  of  these  films  has  been  observed  both  in  situ  and  also  ex  situ. 
Very  recently,  scanning  tunnel  microscope  (STM)  experiments  [10,  11]  have 
shown  phase  transitions  of  a  complicated  nature  occur  in  the  underpotential 
deposition  of  copper  on  gold  (111)  surface.  As  was  hinted  by  the  extended 
X-ray  absorption  fine  structure  (EXAFS)  experiments,  [12],  the  anion,  in 
l liia  case  the  sulfate,  plays  an  important  role  in  the  phase  formation. 


n  For 

The  theory  of  the  UPD  [4,  13]  has  adressed  the  question  of  the  binding 
of  individual  atoms  to  the  metal  surface,  but  has  not  considered  the  coop¬ 
erative  effects  of  fhe  atoms  in  the  adsorbed  film.  In  the  present  work  we  ' ! 
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on  a  triangular  lattice,  which  corresponds  to  the  (111)  face  of  an  fee  crystal, 
or  the  (100)  face  of  an  hep  crystal.  We  propose  a  Pade  approximant  that 
interpolates  the  adsorption  isotherm  between  the  Langmuir  isotherm  when 
the  lateral  interactions  of  the  adsorbate  are  turned  off,  and  the  exact  high 
and  low  temperature  series  for  the  Ising  model  [14]  when  the  interactions 
are  on.  The  coefficients  of  the  Pade  approximant  are  generated  by  a  recur¬ 
rence  relation  to  any  arbitrary  order  in  the  fugacity  z ,  and  therefore  can  be 
generated  for  other  lattices  and  multiadsorption  isotherms.  One  drawback 
of  our  formalism  is  that  it  does  not  satisfy  automatically  the  symmetry  [15] 
relationship  between  the  high  and  low  temperature  series.  These  relations 
will  be  satisafied  if  we  take  a  large  number  of  terms  (perhaps  20  for  the  tri¬ 
angular  lattice),  but  we  can  either  use  a  switching  function  to  force  duality 
on  our  expressions,  or  alternatively  use  a  two  point  Pade  which  will  satisfy 
this  relation.  We  give  a  brief  discussion  of  this  latter  approach,  which  has 
the  disadvantage  of  requiring  the  inversion  of  an  N  by  N  matrix,  where  N  is 
the  order  of  the  approximation.  We  believe  that  the  switching  function  Pade 
approximation  produces  good  results  even  using  as  few  as  8  terms. 

In  previous  work  we  have  proposed  a  simple  model  of  chemisorption  on 
crystal  surfaces  [26,  16,  17],  in  which  the  binding  process  of  individual  atoms 
is  essentially  decoupled  from  the  cooperative  effects  that  take  place  at  the 
surface.  This  model,  which  we  call  the  sticky  site  model  (SSM)  is  a  combi¬ 
nation  of  the  sticky  potential  model,  first  proposed  by  Boltzmann  [18].  but 
solved  and  discussed  in  detail  by  Baxter[19],  and  the  adsorption  site  model 
of  Langmuir  [20].  The  sticky  potential  of  Baxter  was  used  by  Perram  and 
Smith  [21]  to  study  adsorption  onto  a  smooth  interface.  In  that  work  only 
the  density  profile  near  a  smooth  hard  wall  was  discussed,  and  the  properties 
of  the  adsorbed  film  were  not  studied.  The  properties  of  this  adsorbed  film 
is  the  central  issue  of  our  present  work. 

The  sticky  potential  has  the  form 

e-[/J«(r »  =  1  + A6(r- O.  (1) 

where  p  =  1/kT  is  the  usual  Boltzmann  thermal  factor,  u( r)  is  the  inter- 
molecular  potential.  A  is  the  stickiness  parameter,  r  —  (.r,  »/•  ~)  is  the  relative 
position  of  the  center  of  the  molecules,  and  a  is  the  diameter  of  the  molecules. 
The  right  hand  side  term  represents  the  probability  of  two  molecules  being 
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stuck  by  the  potential  u( r):  this  occurs  only  when  the  two  molecules  are  in 
contact,  and  for  this  reason  we  use  the  dirac  delta  function  6(r  —  cr_),  which 
is  zero  when  the  molecules  do  not  touch,  is  infinity  when  they  do,  but  the 
integral  is  normalized  to  one.  The  parameter  A  represents  the  the  probability 
of  adsorption  of  an  individual  molecule  onto  the  surface.  It  has  the  form  of 
a  fugacit  v, 

A  =  e-[i#1-4(r)],  (2) 

where  /j,a  is  the  free  energy  of  the  adsorption  of  an  individual  molecule  onto 
the  surface.  This  quantity  is  obtained  from  a  quantum  mechanical  calcula¬ 
tion. 


The  Langmuir  adsorption  sites  can  be  represented  by  a  collection  of 
sticky  sites  of  the  same  form  as  was  suggested  by  Baxter.  Only  that  now 
we  do  not  have  a  sphere  covered  uniformly  by  a  layer  of  glue,  but  rather 
a  smooth,  hard  surface  with  sticky  points,  which  represent  adsorption  sites 
where  actual  chemical  bonding  takes  place.  In  this  model,  Eq.  1  has  to  be 
changed  to 

e-^“'(r)}  =  1  +  \{R)6(z).  (3) 

with 

A(R)  =  ]T]  A6{R-  riia.1  -  n2a2),  (4) 

ni,n2 

Here  R  =  (,r ,  </)  is  the  position  at  the  electrode  surface,  and  z  the  distance 
to  the  contact  plane,  which  is  at  a  distance  a/ 2  from  the  electrode.  In  4, 
ni,n2  are  natural  numbers,  and  ai,a2  are  lattice  vectors  of  the  adsorption 
sites  which  are  the  sites  of  a  lattice  A  on  the  surface.  The  requirement  of 
point  adsorption  sites  rather  than  extended  regions  around  the  sites  is  not 
essential  to  our  discussion.  It  is  clear  that  phase  transitions  will  occur  even 
in  smooth  surfaces,  simply  because  the  two  dimensional  gas  does  ondergo 
such  phase  transitions.  Less  localized  forms  of  the  adsorption  potential  are 
equivalent  to  the  soft  potential  models  discussed  by  Steele  and  collaborators 
[22],  and  models  of  commensurate-incommensurate  adsorption  such  as  the 
ones  discussed  by  McTague  and  Novaco  [23],  and  by  Shiba  [24]. 

Consider  the  case  of  a  simple  salt  dissolved  in  water,  near  a  metallic 
electrode:  in  the  general  case  there  will  be  three  species  of  the  liquid  phase 
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that  will  be  adsorbed  onto  the  metal  surface:  the  anion,  the  cation  and  the 
solvent.  If  the  metal  is  charged,  then  only  the  ions  of  opposite  charge  will  be 
adsorbed  and  will  compete  with  the  solvent  for  the  adsorption  sites.  If  the 
metal  is  neutral,  then  all  three  species  can  compete  for  the  sites.  This  case 
corresponds  to  a  three  state  Ising  model  in  two  dimensions:  it  leads  to  a  very 
rich  and  complicated  phase  diagram  for  the  adsorbate  and  was  discussed  in 
detail  recently  [25].  When  the  metal  is  charged,  then  the  adsorbate  consists 
of  ions  of  the  same  sign,  which  repel  each  other.  This  situation  corresponds 
to  the  two  dimensional  antiferromagnetic  Ising  model,  which  for  first  nearest 
neighbor  interactions  has  an  order-disorder  phase  transition.  In  the  limiting 
case  of  infinite  strong  repulsion,  it  is  mathematically  isomorphic  to  the  hard 
disk  problem  on  a  planar  lattice,  which  is  known  to  undergo  second  order 
phase  transitions.  Alternatively  the  adions  can  be  neutralized  by  the  elec¬ 
trons  of  the  metal,  and  form  metal  bonds  between  the  surface  atoms.  This 
occur  as  in  the  UPD  :  the  interactions  between  neighbors  are  attractive,  as 
in  the  case  of  the  ferromagnetic  Ising  model  on  a  two  dimensional  lattice, 
and  one  should  expect  first  order  phase  transitions  in  the  adsorbed  film.  We 
should  remark  that  the  existence  of  phase  transitions  is  not  a  consequence  of 
the  existence  of  an  ordered  array  of  adsorption  sites:  they  occur  also  in  the 
smooth  surface  as  a  phase  transition  of  the  two  dimensional  gas. 

In  the  limiting  case  of  the  SSM.  the  sizes  of  the  different  species  play 
a  crucial  role  in  the  possible  ordering  of  the  adlayers  at  the  interface.  It 
will  be  convenient  to  picture  the  ion:;  as  having  a  hard  sphere  core  with  a 
diameter  a ,  and  the  solvent  as  having  a  hard  core  with  diameter  an. 
The  lattice  spacing  of  the  atoms  in  the  metal  surface  is  d,  and  because  it  is 
the  most  stable  surface,  we  will  restrict  ourselves  to  the  (111)  surface  of  the 
fee  crystals,  or  the  (100)  face  of  the  hep  crystal,  that  is,  the  triangular  lattice. 
When  the  adions  are  bigger  than  the  lattice  spacing  and  the  adsorption  sites 
are  extended  rather  than  point  like,  then  the  geometry  of  the  adlayer  is 
different  of  that  of  the  substrate,  and  in  general  incommensurate  with  it. 
This  case  was  studied  in  a  separate  communication  [17]. 

In  this  paper  we  consider  the  case  of  a  triangular  lattice  of  sticky  points. 
The  spacing  between  the  sites  of  the  lattice  is  d.  This  surface  is  in  contact 
with  a  solution.  Only  two  states  of  occupation  are  allowed:  the  sites  are 
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either  occupied  by  ail  ion  or  by  the  solvent,  or  alternatively  in  the  case  of  a 
pure  fluid,  by  a  fluid  particle  or  none.  The  fluid  particles  have  an  exclusion 
diameter  cr,  which  may  or  may  not  be  associated  with  a  hard  core  potential. 
Otherwise  the  interactions  are  arbitrary.  We  assume,  however,  that  the  pair 
correlations  on  the  surface  decay  sufficiently  fast  so  that  we  need  to  take  into 
account  first  neighbor  interactions  only.  There  are  two  possible  situations: 
If  the  adsorbate  diameter  a  is  smaller  than  the  lattice  spacing  d,  then  there 
are  two  possible  phases,  a  dense,  crystalline  one  and  a  dilute  disordered  one. 
There  is  a  first  order  transition  between  them.  If  the  adsorbed  particles 
exclude  all  next  nearest  neighbors,  that  is  when 

d  <  a  <  v/3 d,  (5) 

then  the  problem  is  exactly  analogous  to  the  hard  hexagon  problem  of  Bax¬ 
ter  [15].  In  this  case  there  is  a  second  order  phase  transition  between  an 
ordered  sj 3x^3  phase  and  a  disordered  one.  This  interesting  case  will  not 
be  discussed  here. 


2  Basic  Formalism 

Consider  a  fluid  of  only  one  kind  of  particle  of  diameter  <x  .  near  a 
smooth,  hard  wall  with  sticky  sites.  The  fluid  has  N  particles  and  the  volume 
of  the  system  is  V.  The  Hamiltonian  of  the  system  is 

ff  =  Ho  +  Ht.  (6) 

where  Ho  is  the  Hamiltonian  of  the  system  in  the  absence  of  the  sticky  sites 
on  the  hard  wall,  and  Hs  is  the  sticky  sites  interaction 


M 

H3  =  Y,lli{r),  (~) 

i=i 

where  it3(ri)  is  the  sticky  interaction  of  Eq.  3.  The  canonical  partition 
function  of  this  model  is 
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(8) 


z  =  (1/iV!)  / rfrVJ"»  n[l  +  A(R,)6(;,)]. 

J  1=1 

Expanding  the  product  in  Eq.8  and  integrating  the  Dirac  delta  functions 
we  get,  using  the  single  component  notation  to  avoid  heavy  and  unnecessarily 
complex  equations,  (  with  the  understanding  that  in  the  multi  component 
case  N  is  a  vector  quantity  with  components  A2,  ..-V„  ,  the  necessary 
modifications  of  N!  and  the  integrations  have  to  be  made). 

Z  =  ZoE(X'/nl)  £  rf(R,.R2..R»)  (9) 

n=0  {R,  }CA 

where 

pH i.,r...r.)  =  (Z0(N  -  «)!)“■  J dt.+l. (10) 
=  gi(Tur2..rn)  n^ilr.),  (11) 

i=i 

Here  r2..r„)  is  the  n-body  correlation  function,  and  /9j(r,)  =  p?(-,)  is 

the  singlet  density  of  the  inhomogeneous  smooth  wall  problem.  The  partition 
function  is 

Z„  =  (l/JV!)  y  rfrve-JS».  (12) 

In  the  sticky  sites  model  (SSM),  the  excess  properties  of  the  interface 
depend  only  on  the  correlation  functions  of  the  smooth  interface.  Introducing 
the  potentials  of  mean  force  u’n(Ri.  R2..Rn) 

a:(R„RJ..R.)  =e-J“-IR'R’R-',  (13) 

Combining  these  expressions  we  arrive  at 

Z/Z0=  £  [Api(0)]"/n!e“^u’’,(r,'r2  "rn), 

R.CA 

which  is  the  central  quantity  of  our  work.  The  excess  free  energy  is 
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(14) 


=  (-l/3A)!n(Z/Z0). 


(15) 


where  A  is  the  area  of  the  interface.  We  also  deduce  the  the  fraction  of 
occupied  sites  6 ,  [26] 


A3 X  0Afs 
|A|  <9A 


(1C) 


Eq.  14  shows  that  the  SSM  model  maps  the  adsorption  on  a  flat  sur¬ 
face  onto  a  two  dimensional  lattice  problem  of  a  very  general  kind,  which 
is  in  general  not  amenable  to  analytic  treatment.  It  can  be  simplified  by 
introducing  the  Kirkwood  superposition  approximation 


9l(Ri,R!..R.)  =  n  (1') 

i,j<n 

where 

<,2°(R„R,)=3°(|Rj-RJ|).  (18) 

is  the  pair  distribution  function.  Further  restriction  of  the  interactions  to 
nearest  neighbors  makes  this  problem  equivalent  to  that  of  the  Ising  model. 
The  partition  function  can  be  written  as 

z/z°  =  u,j  t,  =  o.  i.  (i9) 


where  3  =  1/kT  , 


,*  =  -  ln[S2"(rf)],  (20) 

gl(d)  =  gl(  R„R,  )  [/.  j  =  nearest  neighbors]  (21) 

=  ln[Ap?(0)]  (22) 

This  partition  function  can  be  mapped  onto  an  Ising  model  with  spin  vari¬ 
ables  S{  =  ±1  by  means  of  the  transformation 

s,  =  2 U  -  1.  (23) 
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When  we  turn  off  the  interactions  in  the  surface  then  e 
becomes 

z/z°  =  £V"£.\ 

<M 


=  0  and  Eq.  19 


and  using  Eq.  16.  we  get 


Furthermore  since 


Z/Z°  =  (1  +  e',")|A|. 


9  =  A~hi[(l  +  eJ")]. 


e*1  =  ApJ(O), 


we  get  the  Langmuir  isotherm 


l  +  A/>?(0) 

3  Phase  Transitions  at  Electrode  Surfaces 

There  is  a  phase  transition  in  the  ferromagnetic  Ising  model  only  when 
the  external  field  is  zero  [27].  For  the  two  dimensional  lattice  gas  this  means 
that  when  it  <  0 

/(  =  wqf  2  (29) 

or 

V(o)  =  [S2vr,/2-  oo) 

where  q  is  the  number  of  nearest  neighbors  of  the  lattice.  4  for  the  square 
lattice  and  6  for  the  triangular  lattice.  The  exact  spontaneous  magnetization 
of  the  ferromagnetic  Ising  model  on  the  triangular  lattice  [28]  yields  the 
expression 


e  (1/2)  II  ±  ^1  (^2_  + 
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Setting  9  =  (1/2)  tliis  equation  yields  the  condition 

32  =  g\(d)  |rrFf=  3.  (32) 

and  the  value  for  the  critical  sticky  parameter  A  is 

Ap°(0)  |cr„  =  1/27.  (33) 

As  was  mentioned  already  in  the  electrochemical  case,  the  contact  pair 
correlation  function  in  the  bulk  for  ions  of  equal  sign  is  practically  zero,  be¬ 
cause  of  the  Coulomb  repulsion  which  prevents  ions  of  equal  sign  to  approach 
each  other.  In  the  adsorbed  layer  the  interactions  of  these  ions  must  be  at¬ 
tractive,  if  the  formation  of  a  layer  occurs  suddenly,  and  therefore  the  state  of 
chemical  bonding,  or  more  precisely  the  electrovalence,  must  change  during 
the  adsorption  process. 

Condensed  phases  in  the  ad  layers  are  observed  in  electrochemistry.  In 
particular  the  under  potential  deposition  of  some  metals  on  electrodes  occurs 
at  certain  very  well  defined  values  of  the  potential  bias  [8],  For  example,  the 
deposition  of  Cu  on  the  Au(lll)  face  forms  two  phases  according  to  the 
deposition  potential.  These  phases  have  been  observed  ex-situ  [8]  and  in 
situ  [5]  [29]  [30],  [10].  [11].  At  a  lower  potential  a  dilute  ordered  phase  is 
formed.  At  a  higher  potential  a  dense  commensurate  phase  is  formed.  It  is 
clear  front  the  above  considerations  that  in  the  dense  ad  layer  case  the  ions 
must  be  discharged,  because  then  they  would  form  a  metallic  bond,  which 
makes  w  negative,  and  therefore  ferromagnetic.  This  is  supported  by  the 
features  of  the  EXAFS  spectra.  In  the  high  density  phase  the  near  edge 
structure  corresponds  to  that  of  metallic  copper,  wliirh  has  a  characteristic 
double  peak. 

In  the  electrosorption  of  ions,  the  charge  of  the  ions  can  be  neutral¬ 
ized  by  the  electons  in  the  metal  electrode  substrate.  If  this  happens  then 
the  normally  repulsive  effective  interaction  between  equally  charged  ions  can 
become  attractive  because  of  the  formation  of  a  metallic  bond.  The  change 
is  known  as  the  electrosorption  valency,  and  has  been  studied  extensively  by 
Schultze  and  coworkers  [31].  From  the  structure  of  the  Fqs.30  and  Eq.  31 
it  is  clear  that  no  phase  transition  will  occur  if  the  adsorbed  ions  keep  their 
charge,  and  their  repulsive  interaction  as  the  potential  changes.  The  conclu¬ 
sion  is  that  the  adions  ,  in  this  case  Cu.  attract  each  other  in  the  adsorbed 
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layer,  and  therefore  are  chemically  different  in  the  adlayer  than  in  the  bulk 
solution. 

The  contact  probability  parameter  p*( 0)  is  a  function  of  the  electric  poten¬ 
tial:  an  estimate  of  the  contact  can  be  obtained  using  the  Poisson-  Boltzmann 
expression  used  by  the  classic  Gouv-  Chapman  theory  [32] 

Pi(  0)  =  e~s,'tPi ,  (34) 

where  c,  is  the  electrovalence  and  p,  is  the  bulk  density  of  ion  i.  is  the 
adimensional  potential  bias  with  reference  to  the  so  called  potential  of  zero 
charge,  and  is  given  by 

*  =  *fr(0)  -  v'*'].  (35) 

where  e  is  the  elementary  charge,  f3  =  1/kT .  ip( 0)  is  the  potential  at  the 
electrode  surface,  and  ippzc  is  the  potential  of  zero  charge.  Empirically,  recent 
experiments  by  Anderson  [33]  have  shown  that  Eq.  34  works  in  systems  where 
the  assumptions  of  the  Gouy-Chapman  theory  do  not  apply  ,  and  where  we 
know  that  the  discrete  nature  of  the  solvent  will  produce  oscillatory  charge 
profiles  [34].  Yet  contact  theorems  of  the  type 

Wj><  0)  =  +  *T5>,  (36) 

are  valid  for  irregular  surfaces  with  arbitrary  interactions  [12]  in  the  mean 
field  sense.  The  reason  is  that  the  contact  density  p,(0)  given  by  Eq.34,  for  a 
particular  value  of  the  potential  of  zero  charge  U'psc,  is  accurate  is  presumably 
that  the  pair  interactions  cancel  out  at  the  surface,  where  only  the  particle- 
wall  interactions  exist,  and  in  the  mean  field  limit  a  relation  like  Eq.36  is 
satisfied  by  this  ansatz. 

The  fraction  of  occupied  sites  9  Eq.  16  can  be  represented  by  a  series 
[14],  at  both  low  temperature  (liigh  fugacity  z)  and  high  temperature  (low 
fugacity  z)  .  We  consider  the  adsorption  of  a  single  ion  only,  (always  the 
counterion), 


2  =  e~z'*pi\  =  pt(0)A. 

and  in  Domb's  notation 

V  =  «7(*)< 


(37) 

(38) 
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(39) 


a  =  ef3  = 


bV)] 


-l 


9  at  low  fugacity  is  given  by 


1-  0  =  Y^  vyrgr{u ), 

r=  1 


and  at  high  fugacity  is  given  by 


(40) 


#  =  Enrr9r(">-  (4i) 

r=l 

Here  gr{u)  are  polynomials,  obtained  from  graph  theory  [14]  (there  are  not 
to  be  confused  with  the  distribution  functions).  The  first  polynomials  are 


g\{u)  =  u\  (42) 

9i(u)  =  (l/2)(6tt5  -  7u6),  (43) 

g3(u)  =  (l/3)(6u6  +  27 u7  -  90 u8  +  58u9).  (44) 

This  series  must  satisfy  certain  limiting  conditions:  when  y  =  1,  the  critical 
line  should  be  obtained,  for  which  an  exact  expression  is  known.  For  the 
triangular  lattice  Potts  found  [28] 


9  =  (1/2) 


(16u3) 

(1  -  u)3(l  +  3u) / 


(45) 


A  similar  expression  exists  also  for  the  square  lattice.  From  the  high 
and  low  temperature  series  the  correct  critical  exponent  0  =  1/8  can  be 
extracted.  However  the  series  has  poor  convergence  properties.  Furthermore 
Eq.  45  diverges  when  the  lateral  interactions  between  the  adions  is  turned 
off,  that  is,  from  Eq.  38  when  w  =  0  or  u  =  1.  The  observation  is  that  in 
this  limit  both  series  Eqs.  40  and  41  yield  the  series 


e  =  y-1-y~2  +  y-3- 
0  =  l-y  +  y2-y3- 


y  1  _  £ 

1  +  y~ 1  1  +  z 

1  _  2 

ITv  = 


(46) 

(47) 
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The  resummed  series  is  just  Langmuir's  adsorption  isotherm,  which  is  indeed 
the  correct  physical  limit  of  6  obtained  from  Eqs.15,  16  and  19  when  the 
lateral  interactions  are  turned  off.  This  immediately  suggests  the  form  of  the 
Pade  approximant  for  general  2  and  w  =  1/u  =  g2.  For  low  fugacities 


g  _  z  +  z(u'  -  l)P(w.z) 

1  1  4-  2  +  z{iv  -  l)P{a\  z) ' 

and  for  high  fugacities 

a  =  _ I _ 

h  1  +  1/[2U-6]  +  (1/[2U>6])(u'-  l)P(t/',l/[2U)6])’ 


(48) 


(49) 


The  Polynomial  P(u\a)  and  the  details  of  this  Pade  approximation 
are  discussed  in  appendix  A,  where  we  give  explicit  formulas  up  to  order 
8.  As  mentioned  in  the  introduction,  this  approximation  does  not  satisfy 
the  duality  symmetry  relation.  An  alternative  Pade  which  does  satisfy  this 
relation  is  discussed  in  Appendix  B.  Our  approximation  here  will  consist  in 
adding  the  two  expressions  for  the  high  and  low  fugacity  limits  with  a  smooth 
swatch  function,  which  we  have  taken  to  be 

rj{y,s)  -  (1/tt)(^  +  arctan(ys)],  (50) 

where  s  is  .001  for  w  >  3  and  s  =  (3.001  -  w)/b  -when  w  <  3.  This  averaging 
procedure  produces  an  isotherm  that  has  the  correct  shape  for  both  the 
region  above  the  critical  value  of  w  =  3  and  for  the  region  below,  where 
no  sharp  transition  is  expected  to  occur.  In  figure  1  we  show  a  series  of 
isotherms  for  different  values  of  u\  and  we  show,  for  comparison,  the  critical 
line  from  the  exact  solution  of  Potts.  Since  the  calculation  was  done  with 
only  8  terms,  we  think  that  the  agreement  is  satisfactory.  We  could,  with  the 
published  expressions  of  gr(u)  go  up  to  order  20,  but  this  is  certainly  beyond 
the  accuracy  of  our  current  needs.  The  final  expression  for  the  fraction  of 
occupied  sites  is 

0  =  n(y,s)di  +  [1  -  ^(i/.s)]#/,,  (51) 


A  final  observation  is  that  the  mean  field  result  [35]  can  also  be  cast  in 
the  form  of  a  modified  Langmuir  adsorption  isotherm: 


0/  = 


zwq0 

1  4-  zwq9 ' 


(52) 
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We  also  quote  the  results  of  the  two  point  Pade  approximant,  which  is 
discussed  in  Appendix  B.  The  advantage  of  this  approximant  is  that  it  satis¬ 
fies  automatically  the  duality  theorem.  The  disadvantage  is  that  it  requires 
the  inversion  of  an  N  by  N  matrix,  where  N  is  the  order  of  the  approximant. 
We  have  computed  this  Pade  up  to  order  5,  and  found  that  this  is  not  enough 
since  the  adsorption  isotherm  near  the  critical  line  deviates  strongly  from  the 
exact  result  for  the  triangular  lattice.  And  also,  even  the  order  5  Pade  yields 
an  expression  that  takes  about  20  pages  to  write  down.  Above  the  critical 
point  a  [2,  2]  two  point  Pade  approximant  yields 

a  _  (-7  +  6g2  +  ff26)  z  +  6(-l  +g2)g26z 2 
6  (-1  +  <72)  <72622  +  (-1  +  <726)  (1  +  z)  ’ 

The  fraction  of  occupied  sites  as  a  function  of  applied  potential  is  shown 
in  figure  2  for  two  values  of  w  =  2.3,  well  below  the  critical  value  for  film 
formation,  and  for  w  —  3.1,  which  is  above  the  critical  point.  Clearly  a  phase 
transition  occurs  when  y  —  1/ zw 3  =  1.  This  shows  also  in  the  voltamogram 
as  a  sharp  peak.  The  charge  potential  curve  of  the  voltamogram  can  be 
obtained  by  differentiation: 


/(V)  = 


80  dll' 

d^lt' 


(54) 


If  the  scanning  rate  is  constant  and  we  neglect  diffusion  and  double  layer 
effects  [32]  the  derivative  of  0  with  respect  to  the  applied  potential  should 
give  a  repesentation  of  the  voltamogram. 

The  curves  are  shown  in  figure  3,  and  show  a  sharp  peak  for  the  first  case, 
and  a  rather  broad  peak  in  the  case  w  =  2.3  . 


We  believe  that  the  procedure  discussed  in  tills  work  can  be  extended 
to  other  systems,  and  in  particular  to  multicomponent  adsorption,  where  the 
location  of  the  critical  points  is  known,  and  where  sufficient  terms  of  the  high 
and  low  fugacity  series  are  known. 
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4  Appendix  A:  The  Saturation  Pade  Ap- 
proximant 


In  this  appendix  we  discuss  the  Pade  approximant  that  has  the  general  form 


_  A(z) 

1  +  A(z)' 


(55) 


where 


(56) 


i=i 

and  we  have  to  find  the  coefficients  .4,  in  this  equation.  We  remark  first  that 
this  equation  is  the  natural  generalization  of  Langmuir’s  adsorption  isotherm 
to  the  case  of  interacting  films,  since  for  real  and  positive  A(z)  ,  1  >  6  >  0. 
The  low  fugacity  series  Eq.  40  can  be  written  as 


N 

B(z)  =  Y,B,z‘,  (57) 

1=1 

where  the  coeffients  £?,  are  functions  of  m  and  can  be  obtained  directly  from 
the  work  of  Domb  [14]:  The  first  few  are 

Bx  =  1,  (58) 

B2  =  -7  +  6u\  (59) 

Bz  =  58  -  90  m  +  27  m2  +  6 w3,  (60) 

Ba  =  -519  +  1152 w  -  714 w2  +  20u>3  +  48u;4  +  12u>5,  (61) 

B5  =  4856-13980ui+13185tc2-3400u’3-885ud+90ui5+105u;6+30w7,  (62) 


Be  =  -46780  +  165330u>  -  209520m2  +  100842u>3 
-819io4  -  7920tu5  -  1668u;6  +  198m7  +  252 m8  +  84m9,  (63) 


B7  =  460027  -  1926288m  +  3065979m2  -  2152332m3  + 

447090 m4  +  148176m5  -  21049m6 
-18774m7  -  3948m8  +  168m9  +  735m10 

+210mn  +  7m12,  (64) 
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Bg  =  -4593647  +  22229424a  -  42593040 a2  +  39240200 a3  -  15262770a-4 
-469344a5  +  1327156a6  +  246600a7  -  70872a8  -  42504a9 

-14580a10  +  960  a;11  +  1816  a12  +  552a;13  +  48a14(65) 

The  Pade  equation 

+  (66) 

1=1  1=1  j= 1 

yields  a  recurrence  relation  for  the  coefficients  .4,.  We  get 


Ax  =  1, 

(67) 

A2  =  B\  +  B2  —  6(  w  —  1), 

(68) 

A3  =  3  (-1  +  a)  (-15  4-  llu;  +  2 a2)  , 

(69) 

and  so 

on.  After  some  straightforward  algebra  we  find 

Q  z  +  z(w  —  l)P(u\  z) 

1  1  +  z  +  z(w  -  l)P(w,  z)' 

(70) 

where 

N 

P{w,z)  = 

1=1 

(71) 

and 

Pi  =  6, 

(72) 

P2  =  3(— 15  +  11a  +  2  a;2), 

(73) 

Pi  —  2(187  -  266a;  +  46a;2  +  30  a3  +  6u;4), 

(74) 

P*  = 

:  3(-1100  +  2306a;  -  1190a;2  -  156a’3  +  83a;4  +  45a’5  4-  10a’6), 

(75) 

p5  =  6(5047  -  13944a>  +  11897a;2  -  1778a;3  -  1276a;4 

-114a;5  +  99a;6  +  56a7  +  14a8), 

(76) 
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(77) 


p6  =  -285682  4-  978458  a-  -  1161412  a-2  +  459060a-3 
4-67758 a»4  -  41124a-5  -  17876a-6  -  1652a-7  +  1288a-8 

4-952u-9  +  217a10  +  7an, 

p7  =  2751258  -  11237190a  +  17065398a2  -  10735704a;3 
4-11507l6a4  4-  1126218a-5  4-  41376a;6  -  117342a-7 
-43824a8  -  7788a9  4-  3810a10  +  2430a11 

4-600a12  4-  48a13,  (78) 


5  Appendix  B:  The  two  Point  Pade  Approx- 
imant 

Consider  the  lattice  gas  with  first  neighbor  attractions  on  the  triangular 
lattice  .  Following  the  work  of  Domb  [14]  we  have  high  density  and  low 
density  expressions  for  the  excess  free  energy  and  the  fraction  of  occupied 
sites  0  of  the  lattice  gas.  For  low  z  we  have 

0  =  z  4-  (  —  7  4-  65-2)  z2  4-  (58  —  90^2  4-  27g22  4-  65-2^  23....  ( <9) 

,v 

e-A{z)-^Aiz\  (80) 

1=1 

For  high  z 

1—0  =  1/ (2<726)4-(— 7  4-  65-2)  /(2<726)24-(58  -  905-2  4-  275-2 2  4-  65-23)  /  (2<726)  •••• 

(81) 

°r  v 

1  -e  =  B{z)  =  Y,Bilz\  (82) 

»=i 


The  adsorption  isotherm  can  be  constructed  from  both  these  series  by 
a  two  point  Pade  approximant,  which  is  required  to  satisfy  both  Eqs.  5  and 
82  simultaneously.  In  fact,  consider  the  Pade 


a{z) 

14-5(2)' 


(83) 
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where 


a(z )  =  aiZ' ’ 


M*)  = 


which  for  z  — *  0  must  satisfy,  to  order  N 

w  \  a(2) 

V  1  +b(z) 

and  for  z  — ►  oo  must  satisfy  to  order  N 

a(z)/zx 


B{z)  = 


(1  +b{z))/zx' 


(84) 

(85) 

(86) 

(87) 


This  forms  a  system  of  2  N  linear  equations  that  can  be  solved  in  closed 
form:  Defining  the  matrices 

1  Bi  B2 

W  =  -Ax  1  Bx  (88) 

—  A2  —Ai  1 

'1  Bx  B2' 

B  =  0  1  Bx  (89) 

0  0  1 

'-4!  ' 

A  =  A2  (90) 

we  get 

a  =  W~'.A,  (91) 

and 

b  =  B.a.  (92) 

which  are  vectors  of  the  coefficients  of  the  Pade  Eq.  83  a  =  (ax.a2, ..)  and 

6  =  {bx,b2, ...)  . 
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Figure  1:  0  as  a  function  of  the  fugacitv  z  for  w  =  3.5, 3.1, 3.0  (critical 
value)  ,2.95,2.9,2.8,  and  2.6.  The  squares  are  the  spontaneous  transition 
curve  of  Potts. 

Figure  2:  Adsorption  isotherm  as  a  function  of  the  reduced  potential 
'!/,  for  two  values  of  w  =  2.3  and  3.1.  Here  z0  =  0.1. 

Figure  3:  Ideal  voltamogram  with  no  diffusion  effects.  We  show  a 
reduced  intensity  in  units  of  9/ty  for  the  same  parameters  of  Figure  2. 

Figure  4:  The  case  with  w  =  3.1  blown  up  10  times.  This  is  a  crude 
way  to  picture  the  effect  of  diffusion  broadening  on  the  voltamogram. 
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